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Abstract 

"^>^'. In this article, we calculate the mass modifications of the vector and axialvector mesons 

■ D* , B* , Di and Bi in the nuclear matter with the QCD sum rules, and obtain the mass- 
(3 ; shifts SMd- = -71MeV, SMb- = -380 MeV, SMdi = 72MeV, SMb^ = 264 MeV, and the 
^N) ■ scattering lengths a_D* ~ — 1.07fm, as* = — 7.17fm, ani = 1.15 fm and ubi = 5.03 fm for the 

^jr^ D* N, B* N, D\N and B\N interactions, respectively. 

^3 ■ PACS numbers: 12.38.Lg; 14.40.Lb; 14.40.Nd 

, Key Words: Nuclear matter, QCD sum rules 

,__,! 1 Introduction 

The modifications of the hadron properties in the nuclear matter can affect the productions of the 

■ open-charmed mesons and the J/ip in the relativistic heavy ion collisions, the higher charmonium 
Q _ states, such as the ip' , xa, Xc2, etc, are considered as the major source of the J/tp [I]. The charmed 

. mesons can obtain mass augments or reductions in the nuclear matter, if the mass modifications 

I are large enough, the decays of the higher charmonium states to the charmed meson pairs can be 

i facilitated or suppressed remarkably due to the available phase-space, thus the decays to the lowest 

^ ■ state J/ip are greatly modified For example, the higher charmonium states can decay to the 

On . DD pairs instead of decaying to the lowest state J/?A, if the mass reductions of the D and D mesons 

are large enough. On the other hand, the suppression of the J/ip production in the relativistic 
heavy ion collisions is considered as an important signature to identify the possible phase transition 
to the quark-gluon plasma [3]. We should be careful before making definite conclusions. 

The QCD sum rules is a powerful theoretical tool in studying the in-medium hadronic properties 
[1], and has been applied extensively to study the light-flavor hadrons and charmonium states in 
f—i . the nuclear matter [SJ [SI [7] . The works on the heavy mesons and heavy baryons are few, only the 

; D, B, Do, Bo, Ac, Ab, and are studied with the QCD sum rules 111 [TOl El [II US] . The 

. ^ . heavy mesons contain a heavy quark and a light quark, the existence of a light quark in the heavy 

■ mesons leads to large difference between the mass-shifts of the heavy mesons and heavy quarkonia 
^ . in the nuclear matter. The former have large contributions from the light-quark condensates, while 

■ - - ■ the latter are dominated by the gluon condensates [7j [51 [HI [IHl Ej- In this article, we study the 

mass modifications of the vector mesons D*, B* and axialvector mesons Di, Bi in the nuclear 
matter using the QCD sum rules. The present predictions can be confronted with the experimental 
data from the CBM and PANDA collaborations in the future [T41IT5] . 

The article is arranged as follows: we study the mass modifications of the vector and axialvector 
mesons D* , B* , Di and Bi in the nuclear matter with the QCD sum rules in Sec. 2; in Sec. 3, we 
present the numerical results and discussions; and Sec. 4 is reserved for our conclusions. 

2 Mass modifications of the D*, B*, Di and Bi in the nuclear 
matter with QCD sum rules 

We study the mass modifications of the D* and Di mesons in nuclear matter with the two-point 
correlation functions Ilp,^{q), 

Tl^q) = I I d\e''^-^T{j^{x)Jl{0)])p,, (1) 
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where the J^{x) denotes the isospin averaged currents rj^^x) and rj^^^x), 



m,{x) ^ ,t^(,)^ ^»™(^)+^»W(^) , (2) 

which interpolate the vector and axialvector mesons D* and Di, respectively, the q denotes the u 
or d quark. The cq and qc mesons maybe obtain different mass modifications in the nuclear matter, 
just like the if+ and mesons [TB], for example, Hilger et al observe that there exist particle- 
antiparticle mass splittings for the scalar and pseudoscalar mesons [lOlIll] . In this article, we intend 
to study whether or not the decays of the higher charmonium states to the D*D* and DiDi states 
are facilitated in the phase-space, and prefer the average values as the particle-antiparticle mass 
splittings cannot modify the total mass of the particle-antiparticle pair. 
At the low nuclear density, the in-medium condensates (0)p„, 

(0)p„ - <^) + 2^<'^)^' 

based on the Fermi gas model, where the (O) and (O) at denote the vacuum condensates and 
nuclear matter induced condensates, respectively, the pat is the density of the nuclear matter [5]. 
Accordingly, the correlation functions 11^^(5) can be divided into a vacuum part 11^^(9) and a 
static one-nucleon part (q) in the Fermi gas approximation for the nuclear matter [51 [5] , 

np.(g) = n°,(g) + 2^T^^,(g), (4) 

where 

T^.i^,q) = *|dW-^(A^(p)|r{jp(a:)jt(0)}|iV(p)), (5) 

the \N{p)) denotes the isospin and spin averaged static nucleon state with the four-momentum 
p = {Mn,0), and normalized as {N{p)\N{p')) = {2TTf2po5^{p-p') [S]. The T^^{q) happen to be 
the current-nucleon forward scattering amplitudes. We can decompose the correlation functions 
T^A^,q) as 

Tj:(c.,q) = T^(^u;,q)(^~g,^+^-I^^+U%{u;,qfj^, (6) 

according to Lorentz covariance, where the Tjv(w, q ) denotes the contributions from the vector and 
axialvector mesons, and the T^{uj,q) denotes the contributions from the scalar and pseudoscalar 
mesons. The interpolating currents '7^(2;) and 775^ have non- vanishing couplings with the scalar and 
pseudoscalar mesons £>o and!?, respectively, i.e. (0|77p(0)|-Do+^o) = fDoQ^i and {0\r]5f^{0)\D+D) = 
ifoqfi- We can exclude the contaminations by choosing the tensor structure —Qf^u + and only 
take into account the vector and axialvector mesons through the definitions, 

{0\7^^mD* + D*) = fD'MD-e^, 

(0|775p(0)|i5i+^i> - fD.Mo.e^, (7) 



9^ 



with summations of the polarization vectors £^('^7 q) — — .g^jy 

In the limit of the 3- vector 7 — >■ 0, the correlation functions Tn{uj, q ) can be related to the D*N 
and DiN scattering T-matrixes, i.e. Td'n{M^d* ,Q) = 87r(Afjv -I- MD*)aD* and Tdxn{Mdi,Q) = 
87r(Afjv + A/£iJa£ii, where the ao' and a^i are the D*N and DiN scattering lengths, respectively. 
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Near the pole positions of the D* and Di mesons, the phenomenological spectral densities p{uj, 0) 
can be parameterized with three unknown parameters a, h and c 0, 



p(w,0) = 



Td'/Din{^^:0) 



le] 



JJ 



(8) 



(9) 



the terms denoted by • • • represent the continuum contributions. The first term denotes the 
double-pole term, and corresponds to the on-shell (i.e. = effects of the T-matrixes, 



and related with the mass-shifts of the D* and D\ mesons through the relation 



Pn 



■a; 



(10) 



(11) 



the second term denotes the single-pole term, and corresponds to the off-shell (i.e. 7^ ^d*/Dx) 
effects of the T-matrixes; and the third term denotes the continuum term or the remaining effects, 
where the so is the continuum threshold. 

In the limit cj — > 0, the Tisi{uj, 0) is equivalent to the Born term T^Jy'^^^(w, 0). We take into 
account the Born term at the phenomenological side. 



So - 



(12) 



with the constraint 



C 
So 



— = 0. 



(13) 



The contributions from the intermediate spin-| charmed baryons are zero in the soft-limit 
— > ^17,, and we only take into account the intermediate spin-i charmed baryons in calculating 
the Born terms, and parameterize the hadronic matrix elements as 



{A,/^,ip~q)\D*i-q)N{p)) - UajsAp-i) 
{AJj:,{p-q)\D,{-q)Nip)) = UAj^Jp~q) 



M, 



N 



Ac/Sc 



gA,/s, ^+ z a Pe^^qp 

MN + JWA JSe 



15Un{p) , 
(14) 



where the Un and Ua^/t.^ are the Dirac spinors of the nucleon and the charmed baryons Ac/Sc, re- 
spectively; the gAjs^D'N, ffAe/Se-Diiv, 9ajt.,D'N ^^'^ sL/s.DiAt ^he strong couphng constants 
in the vertexes. In the limit (7^ — > 0, the strong coupling constants ^-^ j-^.^^ 
no contributions. 

We draw the Feynman diagrams, calculate the Born terms and obtain the results 



and glj^^D.N have 



pBorn , 



u;,0) = 



2flMlMM{MH - MN)glD^N 



(15) 
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where the H means either A+, E+, or 1]°. The masses Ma, = 2.286 GeV and Afs, = 

2.454 GeV from the Particle Data Group [18], we can take Mh ~ 2.4 GeV as the average value. 
On the other hand, there are no inelastic channels for the D*N and DiA^ interactions in the case 
of the charmed mesons cq. 

The scattering state D*N can translate to the scattering states D*N, ttEc, f]^c, DN, ttAc, 
pSc, pAc, etc, we can take into account the infinite series of the intermediate baryon-meson loops 
with the Bethe-Salpeter equation to obtain the full D*N — >■ D*N scattering amplitude, and the 
higher resonances, such as the Ac(2595), Sc(2800), etc, appear as dynamically generated baryon 
states [H]. We can saturate the full D*N -> D*N scattering amplitude with the tree-level Feyn- 
man diagrams of the exchanges of the higher-resonances including the negative spin-^ charmed 

baryons Ac(2595) and ^0(5 ) have the average value AIh' ~ 2.7 GeV [T51 [50] and other excited 
spin-i charmed baryons. The translations of the scattering state D*N to the ground states Ac 
and Sc are greatly facilitated in the phase-space, furthermore, the mass-shift SMd* does not sen- 
sitive to contributions of the ground states Ac and Sc, see Table 1; neglecting the loop-effects (or 
higher-resonance contributions) cannot impair the prediction power remarkably, in other words, 
the contaminations from the intermediate hadron states are small. As the scattering amplitudes 
contain both elastic and nonelastic parts, we use the scattering lengths and Born terms to pa- 
rameterize the elastic and nonelastic scattering amplitudes respectively. The contributions of the 
scattering states (or continuum states) ttSc, TyAc, DN, ttAc, pSc, pAc, etc, are embodied in the 
Born terms. 

We carry out the operator product expansion to the condensates up to dimension-5 at the large 
space-like region in the nuclear matter, and obtain the analytical expressions of the correlation 
functions at the level of quark-gluon degree's of freedom. 



(16) 



where the Cn{qQ,q) are the Wilson coefficients, the in-medium condensates {On)pM = {On) + 
2^^(Cri)Af at the low nuclear density, the {On) and {On)N denote the vacuum condensates and 
nuclear matter induced condensates, respectively. One can consult Refs.[SJ |B] for the technical 
details in the operator product expansion. Then we collect the terms proportional to pN (or the 
nuclear matter induced condensates) , take the quark- hadron duality. 



qt^qu 



Y,Cn{^.q){0„ 



>N 



(17) 



set cj^ = 9^, and perform the Borcl transform with respect to the variable — — w^, finally obtain 
the following two QCD sum rules: 

>J- So 



1 



M4. 



e 



bie 



So 
Ml. 



N 



1 



{Mh + Mn)"^ ~ Ml. 



1 

M2 



1 



{Mh + Mn)^ - Ml. 
(Mjj+Mjv)2 1 f mc{qq)N 



{Mh + Mn)^ 



Ml. 



2{qUDoq)N , ml{qUDoq)N , mc{qgs(jGq)N , 8mc{qiDoiDoq) n ml{qiDoiDoq) n 



1 ,a,GG, 



24 



Af2 



N I dx\l 



2Af2 



3M2 



48A/2 



)N 



Af4 

~2 ™c 



(18) 
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So -^1 , 2flMlMN{MH-MN)gln,N 



Ml 



■e > + 



1 



{Mh - MnY - Ml 
(Mfj-Mj^f I (nic{qq)N 



(Mh - MmY -Ml j ~ t 2 

2(qkDQq)N , m'^(qUDoq)N mc{qgs(jGq) n SmdqiDQiDoq) n ml{qiDoiDoq) n 



3Af2 



M4 



1 ,asGG, 



24' TT 



TV / dx\l 



ml \ 1 ,a,GG, 



1 ~ X m'; 



(19) 



where = — . 

a; 

Differentiate above equation with respect to -py , then ehminate the parameter b, we can obtain 
the QCD sum rules for the parameter a. With the simple replacements rric — > mf,, Mu* — >■ Mb*, 
Moi — > Mbi, Ac — > Ab and Sc — > Sb, we can obtain the corresponding QCD sum rules for the 
mass-shifts of the B* and Bi mesons in the nuclear matter, where we take the approximation 



M 



H 



5.7 GeV 



The present approach was introduced by Koike and Hayashigaki to study the spin-isospin 
averaged meson- nucleon scattering lengths and the relevant mass-shifts for the p, u), 4> mesons in 
the nuclear matter [21j . The heavy mesons contain a heavy-quark and a light-quark. The existence 
of a heavy quark in the heavy mesons results in much difference between the in-medium properties 
of the heavy mesons and light mesons. The heavy quark interacts with the nuclear matter through 
the exchange of the intermediate gluons and the modifications of the gluon condensates in the 
nuclear matter are mild, while the modifications of the quark condensates in the nuclear matter 
are rather large. The relation in Eq.(3) is hold for the gluon condensates. We expect that the 
convergent behaviors of the heavy-light type interpolating currents are better than that of the 
light- light type interpolating currents, if the correlation functions n^,y(g) are expanded in terms of 
the external parameter pjv- The approach developed for the light mesons still works for the heavy 
mesons. 

In Ref. 22 , Klingl, Kaiser and Weise use an effective Lagrangian which combines chiral SU{3) 
dynamics with vector meson dominance to calculate the forward vector-meson-nucleon scattering 
amplitudes, and take them as input parameters in the hadronic side of the QCD sum rules in nuclear 
matter, and observe a remarkable degree of consistency with the operator product expansion at the 
quark level. In Ref. [23], Leupold and Mosel study the electromagnetic current-current correlation 
functions in the nuclear matter, and expand the QCD sum rules in terms of the finite squared 
three-momentum g^, and observe that the QCD sum rule can provide an interesting and non- 
trivial consistency check for the hadronic models, but cannot rule out the hadronic models which 
predict a different behavior of the vector mesons with different (7^. In this article, the hadronic 
side of the QCD sum rules of the order 0{pn) consists of both elastic and nonelastic scattering 
amplitudes, which are parameterized by the scattering lengths and the Born terms respectively, 
can lead to stable QCD sum rules with the suitable Borel parameters in a finite range, i.e. the 
predictions make sense. 



3 Numerical results and discussions 

In calculations, we have assumed that the linear density approximation is valid at the low nuclear 
density, (C')p„ = (0|O|0) + ^{N\0\N) = (O) + ^{0)n for a general condensate (0)p„ 
in the nuclear matter. The input parameters are taken as {qq)N — >» Trwd C^^n), { "'^^ )n = 
(-0.65 ± 0.15) GeV(2MAr), {qUDoq)N = (0.18 ± 0.01) GeV(2A/7v), (qgscrGq) = 3.0 GeV2(2AfAr), 
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{qiDoiDoq)N + \{qgs(jGq)N = 0.3 GeV^(2Afjv), m„ + 12MeV, (Jn = (45± 10)MeV, Afjv = 

0.94 GeV, pN = (0.11 GeV)^ [5], = (1.35 ± 0.1) GeV and m,, = (4.7 ± 0.1) GeV at the energy 
scale ^ = 1 GeV. 

For the well established vector mesons D* and B* , we take the values from the Particle Data 
Group, Md* — 2.01 GeV and Mb* — 5.325 GeV [18], the decay constants fo* and fs' are deter- 
mined by the QCD sum rules, fo- = 0.270 GeV and Jb* = 0.195 GeV [M], where the threshold 
parameters are taken as s^,. — (5 — 7) GeV^ and s%, — (33 — 37) GeV^ [24l |25], here we have 
neglected the uncertainties of the decay constants. We can take the threshold parameters as 
s%. = 6.5 GeV^ and s%. = 35 GeV^ to reproduce the values Mo- = 2.01 GeV, Mb- = 5.325 GeV, 
/d* = 0.270 GeV and Jb* = 0.195 GeV approximately for the QCD sum rules in the vacuum. The 
mass of the axialvector meson 01(2430) is Md^ = (2427 ± 26 ± 25) MeV from the Particle Data 
Group [18], and the axialvector meson Bi has not been observed yet. We calculate the hadronic 
parameters of the axialvector mesons Di and Bi using the QCD sum rules in the vacuum, and 
obtain the values Md, = 2.42 GeV, Mb, = 5.75 GeV, fo, = 0.305 GeV and /s, = 0.255 GeV with 
the threshold parameters = 8.5 GeV^ and = 39GeV^. The value Md, = 2.42 GeV repro- 
duces the experimental data AId, = (2427 ± 26 ± 25) MeV well [18 . The prediction of the mass 
Mbi satisfies the relation Mb, — Mb- ~ Md, — Md'- In this article, the threshold parameters 
are taken as s?,. = (6.5 ± 0.5) GeV^ s\ = (8.5 ± 0.5) GeV^ s%, = (35 ± 1) GeV^ and = 
(39 ± 1) GeV^, respectively, which satisfy the relations s^j. Bi — {Md',Di,Bi + 0.4 ^ 0.6 GeV)^ 
and s^. = [Mb' + 0.5 ^ 0.7 GeV)^. In general, the energy gap between the ground state and the 
first radial excited state is about 0.5 GeV. For the explicit expressions of the QCD sum rules in the 
vacuum derived from the correlation functions n"j^(q), one can consult Ref. [26j and the references 
therein. 

The value of the strong coupling constant goNA^ is Qk^dn — 6.74 from the QCD sum rules 
[77] . while the average value of the strong coupling constants qk^dn and g^^oN from the light-cone 
QCD sum rules is g^^eJ^+g^cD" = 6.775 [21], those values are consistent with each other. The 
average value of the strong coupling constants qa^d-n and g^^wN from the light-cone QCD sum 
rules is ^'^°°*"^^^'=°''" = 3.86 [55]. In this article, we take the approximation g\^D' n ~ g^cD* n ~ 

5Ac-DiAf ~ gSaDiN ~ gAbB'N ~ gSbB* N ~ gA^BiN ~ gs^BiN ~ 3.86. 

In Fig.l, we plot the mass-shifts 5M versus the Borel parameter at large intervals. From 
the figure, we can see that the values of the mass-shifts are rather stable with variations of the 
Borel parameter at the intervals M^ = (4.5 - 5.4) GeV^ (6.5 - 7.6) GeV^ (22 - 24) GeV^ and 
(34 — 37) GeV^ for the D*, Di, B* and Bi mesons, respectively; in other words, the uncertainties 
originate from the Borel parameter M'^ are less than 1%. The main contributions come from the 
terms zLmc{qq)N and zLnib{qq) n , see Eqs. (18-19) and Fig. 2, the operator product expansion is well 
convergent. The spectral densities at the level of the quark-gluon degrees of freedom consist of the 
medium-induced condensates, and have the form AiS{s — m'^^i^) + A25{s — fn^^f^), where the Ai and 
A2 denote the coefficients, we carry out the integrals. 



So 



ds 



A^Sis ~ mi/,) + A^Sis ~ mi/,) e"^ , (20) 



to obtain the right side of the QCD sum rules in Eqs. (18-19), there are no perturbative terms to 
approximate the continuum contributions at the regions s > sq. At the phenomenological side, the 
exponential factors 

e-i^ ^ g-(1.20-1.44)^ g-(1.12-1.31)^ g-(1.46-1.59)^ ^-(1.05^1.15)^ (21) 

at the intervals = (4.5-5.4) GeV^ (6.5-7.6) GeV^ (22-24) GeV^ and (34-37) GeV^ for the 
D*, Di, B* , Bi mesons respectively, where we take the central values of the threshold parameters, 
the corresponding exponential factors of the ground states are 

e~l^ = g-(0.75-0.90)^ g-(0.77-0.90)^ ^-(1.18-1.29)^ ^-(0.89-0.97)^ ('22) 
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where the m stands for the D*, Di, B*, Bi mesons respectively; the continuum contributions are 
suppressed more efficiently. Furthermore, we expect that the couplings of an special interpolating 
current to the excited states are more weak than that to the ground state mesons. For example, 
the decay constants of the pseudoscalar mesons 7r(140) and 7r(1800) have the hierarchy: /^(isoo) 
/ir(i40) from the Dyson-Schwinger equation [29], the lattice QCD [30], the QCD sum rules [31], etc, 
or from the experimental data [32| . 

We can take the Borel windows as = (4.5 - 5.4) GeV^ (6.5 - 7.6) GeV^ (22 - 24) GeV^ and 
(34 — 37) GeV^ for the D*, Di, B* and Bi mesons, respectively, and obtain mass-shifts SAId* = 
-7lj:23MeV, SMb* = -380;^i MeV, SMd, = 72tllMeY, SMb^ = 264;^^MeV, respectively; 
and the scattering lengths a^. = -1.07lo;34 fm, as* = -7.17^ J;^! fm, = 1.15tQ;32fm and 
aBi — 5.03t^ 3^fm for the D*N, B*N, DiN and BiN interactions, respectively. For the technical 
details in analyzing the uncertainties, one can consult Ref. |33j . 

In Fig. 3, we plot the mass-shifts dA4 versus the Borel parameter and the strong coupling 
constants g^. From the figure, we can see that the mass-shifts decrease (increase) monotonously 
with increase of the squared strong coupling constants for the vector mesons D* and B* (ax- 
ialvector mesons Di and Bi) in the Borel windows. The precise values of the mass-shifts and 
scattering lengths are presented in Table 1. 

Although the present QGD sum rules are not stable with variations of the Borel parameters, 
the uncertainties originate from the Borel parameters are less than 1% in the Borel windows, i.e. 
we choose suitable platforms to avoid large uncertainties. On the other hand, we can take moments 
of the correlation functions and derive QCD sum rules to study the mass-shifts as in Ref. 1 10), the 
present QCD sum rules are equal to that derived from the first two moments of the correlation 
functions, i.e. we take the weight functions to be 1 and s. 

In the present work and Refs. [H |9] , the correlation functions are divided into a vacuum part 
and a static one-nucleon part, and the nuclear matter induced effects are extracted explicitly, while 
in Refs.[ini E], the pole terms of the hadronic spectral densities are parameterized as ^'"'^("^Q) = 
F+6iLU - M+) - F^d{Lu + M_), where M± = M ± AM and F± ^ F ± AF, and QCD sum 
rules for the mass center M and the mass splitting AM are obtained. For the pseudoscalar D, 
D mesons, Hayashigaki obtains the mass-shift SAId = — 50MeV [8], while Hilger, Thomas and 
Kampfer obtain the mass-shift SMd = -|-45MeV [10]. For scalar Do, Do mesons, the mass-shift 
SAIdo = M — Mdq < obtained by Hilger and Kampfer [11] differs from the result SAIdq = 
-1-69 MeV obtained by Wang and Huang [9]. In Ref. [12], Hilger, Kampfer and Leupold study the 
chiral partners of charmed mesons in the nuclear matter, and focus on the differences between the 
pseudoscalar and scalar as well as vector and axialvector D mesons and derive the corresponding 
Weinberg type sum rules, while the mass-shifts are not presented. 

In the limit iriq 0, the quark condensate (qq) serves as the order parameter and indicates 
that the chiral symmetry is broken. The quark condensate undergoes reduction in the nuclear 
matter, (g(?)p„ = (qq) + 2a"jv (^'?)^' chiral symmetry is partially restored, for example, the 
in-medium nucleon mass Af^ can be approximated as Af^ — —jj^{qq) ^^J^d the mass reduction 
is rather large; on the other hand, there appear new medium-induced condensates, for example, 
the {qiDoiDoq) N , {qgs<yGq) n , etc, which also break the chiral symmetry. In the present case, the 
medium- induced condensates are associated with the large heavy quark masses mq, ttiq, TOq or 
ttiq, the net effects do not always warrant that the chiral symmetry is monotonously restored with 
the increase of the density of the nuclear matter. The light vector current q^y^q and axialvector 
current qjp,j5q are invariant under the chiral transformation q — > e'"^^q, however, the heavy vector 
current Qj^q and axialvector current 07^75? are mixed with each other under the transformation, 
the heavy quark currents Qjf^q and (57^759 are not conserved in the limit niq 0, it is better 
to take the doublets {D*,Di), {B*,Bi) as the parity-doublets rather than the chiral-doublets. If 
we take into account the flavor SU{?>) symmetry of the light quarks, the chiral SU{3)l x SU{3)ii 
transformations require that the ground states (1)*°, D*~) and (B*+, 5*°) have their 
chiral partners {D^, D^ , DJ^) and (Bj*", 5°, -B^i), respectively, those parity-doublets are chiral- 
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doublets. When the density of the nuclear matter is large enough, the order parameter {qq)p^ — > 0, 
the chiral symmetry is restored, the Fermi gas approximation for the nuclear matter does not 
survive, there are free of the non-perturbative contributions from the condensates, and the parity- 
doublets (or chiral-doublets) maybe have degenerated masses approximately. In the present case, 
we study the parity-doublets (or chiral-doublets) in the low nuclear density, the mass breaking 
effects of the parity-doublets (or chiral-doublets) maybe even larger. 

The axialvector current u{x)j^j5c{x) interpolates the axialvector meson £'i(2430) has non- 
vanishing coupling with the scattering state D*tt, in the soft tt hmit, the coupling constant Xd*tt 
can be estimated as 

{0\u{xh^roc{x)\D*7r) = -^{0\[Q5,uixh^,j5cix)]\D*) 

= -^{0\d{xMx)\D*)=-'-^^^^^=zXD'.e^, (23) 

= fD,Mn,e^ ^ \ {D,\D*^), (24) 

where the axial-charge Q5 = J cfiyd^ {y)'^5u{y), and the and are the polarization vectors of 
the vector and axialvector mesons D* and Di, respectively |33], the formula in Eq.(24) survives 
beyond the soft tt limit. The coupling constant Xd'tt is a large quantity and cannot be neglected. 
The rescatterings 

D*TT D*1T, 

D*TT D*T:,D*K,D*rj ^ D*t: , 

D*TT D*T:,DlK,D*rj ^ D*T:,D*K,D*rj ^ D*n , 

D*n D*Tr,DlK,D*r]^ D*TT,D;K,D*r]-^ D*Tr,DlK,D*f]-> D*Tr, (25) 

• • • also have contributions to the hadronic spectral densities. In the heavy meson chiral unitary 
approach, we can use the Bethe-Salpeter equation to perform the summation of the infinite series 
of the intermediate meson-loops (such as the D*Tr, D^K, D*ri) to obtain the full £)*7r — >■ D*7r scat- 
tering amplitude, and generate the axialvector meson Z)i(2430) dynamically [151 [33]. If we saturate 
the full £'*7r D*7t scattering amplitude with the exchanges of the intermediate axialvector meson 
Z?i(2430), the D*tt rescattering effects lead to the renormalization — — —. > — — ^ ^, , — 

in the hadronic representation of the correlation functions, where the n(p) denotes the renormal- 
ized self-energy of the intermediate D*tt loops, and contributes a finite imaginary part to modify 
the dispersion relation. In fact, the contributions of the intermediate meson-loops are very large, 

o 

we have to take the mass Mjj-^ as the bare mass M Di to absorb the real part of the un-renormalized 
self-energy to reproduce the physical mass, the net effects are embodied in the finite imaginary 
part. We can take into account those meson- loops effectively by taking the following replacement 
for the hadronic spectral density, 

6{s-MlJ ^ , (26) 

here we neglect the complicated renormalization procedure for simplicity |36' . In Ref . [36' , we 
observe that a width about (or less than) 400 MeV cannot change the prediction significantly, 
the i5 function approximation for the spectral densities still survives. In the present case, F^i^ = 
3841^°'' ± 75 MeV from the Particle Data Group [T5], the contaminations from the intermediate 
state D*Tr are expected to be small. Analogical discussions can be applied to the contaminations 
from the intermediate state B*tt. 

The negative scattering lengths ajj* — — 1.07fm and qb' = — 7.17fm indicate that the D*N 
and B*N interactions are attractive, it is possible to form the D*N and B*N bound states; while 
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9' 





10 


20 


30 


40 


50 


SMd' (MeV) 


-75 


-72 


-70 


-67 


-65 


-62 


SMb' (MeV) 


-382 


-381 


-380 


-380 


-379 


-378 


5Md, (MeV) 


70 


71 


73 


74 


76 


78 


6MB^ (MeV) 


262 


263 


264 


265 


266 


267 


ao' (fm) 


-1.13 


-1.09 


-1.05 


-1.02 


-0.98 


-0.94 


as* (fm) 


-7.20 


-7.18 


-7.17 


-7.15 


-7.14 


-7.13 


flDi (fm) 


1.11 


1.14 


1.16 


1.19 


1.21 


1.24 


asi (fm) 


5.00 


5.02 


5.04 


5.05 


5.07 


5.09 



Table 1: The mass-shifts SM and the scattering lengths a versus the strong coupling constants 
9'. 

the positive scattering lengths — 1.15 fm and as^ — 5.03 fm indicate that the DiN and BiN 
interactions are repulsive, it is difficult to form the DiN and BiN bound states. According to the 
observations of the present work and Refs.[8l[9], we can draw the conclusion tentatively that the 
negative parity heavy mesons decrease their masses in the nuclear matter while the positive parity 
heavy mesons increase their masses. The decays of the high charmonium states to the negative 
parity charmed meson pairs DD and D*D* are facilitated, while the decays to the positive parity 
charmed mesons pairs DqDq and DiDi are suppressed. The J production can obtain additional 
suppressions due to mass modifications of the negative parity charmed mesons D and D* in the 
nuclear matter. 

4 Conclusion 

In this article, we calculate the mass-shifts of the vector and axialvector mesons D* , B* , Di and 
Bi in the nuclear matter using the QCD sum rules. We take the linear approximation at the 
low density of the nuclear matter, and extract the mass-shifts and scattering lengths explicitly, 
SMd' ^ -71 MeV, SMb' = -380 MeV, SMd, = 72MeV, SMb^ = 264MeV, ao- = -1.07 fm, 
aB' — — 7.17fm, au^ = 1.15fm and a^i = 5.03fm. Our numerical results indicate that the 
D*N and B*N interactions are attractive while the DiN and BiN interactions are repulsive; it is 
possible (difficult) to form the D*N and B*N {DiN and BiN) bound states. The J/^ production 
can obtain additional (no additional) suppression due to mass modification of the vector meson 
D* (axialvector meson Di) in the nuclear matter. The present predictions can be confronted with 
experimental data in the future. 
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